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2 In this article we investigate the possibility of generating piezoelectric or- 

Ch bital polarization in graphene-like systems which are deformed periodically. 

We start with discrete two-band models which depend on control parameters; 

psj in this setting, time-dependent model hamiltonians are described by loops 

>>■ in parameter space. Then, the gap structure at a given Fermi energy gener- 

^P ates a non-trivial topology on parameter space which then leads to possibly 

.^ non-trivial polarizations. More precisely, we show the polarization, as given 

t^ by the King-Smith-Vanderbilt formula, depends only on the homotopy class 

.^ of the loop; hence, a necessary condition for non-trivial piezo effects is that 

^^ the fundamental group of the gapped parameter space must not be trivial. 

Cn The use of the framework of non-commutative geometry implies our results 

. . extend to systems with weak disorder We then apply this analysis to the uni- 

^ axial strain model for graphene which includes nearest-neighbor hopping and 

a stagger potential, and show that it supports non-trivial piezo effects; this is 

in agreement with recent physics literature. 
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1 Introduction and main results 

Piezoelectric materials are crystalline solids which become macroscopically charged when 
subjected to mechanical strain. One material that has recently moved into the limelight of 
piezoelectric physics due to the theoretical work [ ,0R12J is graphene, and an experimental 
realization of these ideas would open up a lot of possibilities in the engineering of new 



piezoelectric devices. Much of graphene's peculiar properties [CGP+09| stem from the 



conical intersections of valence and conduction band right at the Fermi energy. But the 
reason why it is an interesting material for piezoelectric devices is its unique mechanical 
robustness, allowing elastic deformations of up to 20% (as opposed to ^ 0.1 % for normal 
materials) [ |LML07t piVK+OSl |KZJ+09| . 

To understand the link between graphene's band structure and piezoelectric properties, 
one needs a microscopic description of the piezoelectric effect. Such a description had 
eluded theoretical physicists until the mid-1970s when Martin [ Mar74 ] noticed that pre- 
vious definitions of polarization in terms of microscopic dipole moments were incomplete. 



It took another 20 years until Resta l ]Res92[ and King-Smith and Vanderbilt [ ]KV93[ de- 
rived a formula for polarization from linear response theory. They recognized the crucial 
role of the adiabatic Berry phase [ Ber84t IXCN IO*] and linked the difference in charge, 
the polarization AV = [AVi,. ..,AVi), accumulated during a deformation in the time 
interval [0, T] to 



AVj :=i I dtrfPM [d,P[t), V^P(0]). 
Jo 



(1.1) 



Here, T denotes the trace per unit volume, P[t)= 1(-oo,e) (^CO) is the projection onto all 
states below the Fermi energy E and H(t) is the hamiltonian of the system. This equation 
is structurally identical to that for computing other Chern numbers such as those for the 



quantum Hall effect [ TKN-l-82 BES94 1 . One crucial ingredient for the piezoelectric effect 



to occur is the absence of conducting states around the Fermi energy, e. g. materials with 
a spectral gap are good candidates. 

A mathematical justification of dl.lD , also called King- Smith-Vanderbilt formula, has 
first been achieved by Panati, Sparber and Teufel [ PST09 1 for the (commutative) case 



of continuous Schrodinger operators. In a later work, Schulz-Baldes and Teufel [ST13| 
used the language of non-commutative geometry to establish dl.lD for dirty lattice systems, 
i. e. discrete operators which include the effects of random impurities. Both works also 
explore the topological nature of AV in the case of periodic deformations where AV is 
quantized in appropriate unitQ what is missing, however, are criteria that tell us which 
periodic deformations lead to non-trivial polarization AV ^ 0. 

The main focus of this paper is the study of the piezoelectric effect for graphene sub- 
jected to periodic deformations. We will study the simplest kind of tight-binding model, 
the so-called uniaxial strain model. It includes only nearest-neighbor interactions and a 
stagger potential, and will be explained in more detail below. Our investigation has led us 
to the following three questions: 

(Ql) What is the topological origin of non-trivial polarizations? 



^The units for polarization are charge density x distance, and a more careful consideration after restoring 
physical units yields AV = p^r 2],=i ^'Pj Tj where e is the electron charge, |V| the volume of the Wigner- 



Seitz cell and the Yj are a basis for the lattice r |KV93j equation (13)] 
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(Q2) Are there sufficiently general models applicable to graphene for which AV ^ 0? 

(Q3) Is AV stable under perturbations? 

Because our ideas can in principle be applied to any parameter-dependent system, we will 
formulate the first part of this work in more generality. 

1.1 The topology of the parameter space 

For the models we study, the hamiltonian H[q) is described by a set of control parameters 
q = (qi,...,q^,)eR'^, i. e. 

H:Q — *A 

is a continuous (or even more regular) function that takes values in the selfadjoint ele- 
ments of some algebra of bounded operators A. We will always assume that Q is a subset 
of K^. These parameters q model the influence of external effects on the quantum system; 
In our example, the q^ could be hopping parameters. 

Any choice of Fermi energy E singles out configurations denoted with Q^ made up of 
those values of q e Q for which 

(i) E lies in a spectral gap oiH[q), E ^ cr(H(q)), and 

(ii) there are states below E, i. e. P(q) := l(-oo,E)(W(q)) 7^ 0. 

We shall refer to Q^ as the space of gapped configurations at E. It is the topology of Qg 
which determines whether AV = or not; more precisely, the fundamental group ttjCQ^) 
[ |Hat02 ] can provide a classification for the piezoelectric effects for a given model system 
at a given Fermi energy. The idea is as follows: To each physical deformation which does 
not close the gap at E, we can associate a loop in parameter space 17 : §^ ^ Q^ so that the 
time-dependent, T -periodic hamiltonian 

H^(t) = H(i7(2^Vr)). (1.2) 

can be expressed in terms of the loop 17 and the model hamiltonian H. The fact that the 
deformation should be continuous implies that 17 is continuous. Then H^ in turn defines 
a time-dependent Fermi projection 

Pr,(0=l(-oo.£)(H,(0) (1.3) 

which is then plugged into equation ( |1.1| ) to obtain AV[r]) for each loop 17. Overall, this 
procedure yields a map 
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from the space of loops in Qg. 

That AV[r]) is a topological quantity is reflected in the fact that the value depends only 
on the equivalence class [17] e tiiCQ^): if 17 and 17' are homotopic, then also the corre- 
sponding Fermi projections P^ and P^/ can be continuously deformed into one another 



(Proposition 3.21. Thus, the in variance of ( |1.1D under homotopies yields 



Theorem 1.1 (Homotopy-invariance of AP) Under the technical conditions enumerated 



in Theorem 3. 7 the map 

AV,:n,[QE)^Z'^, [17] - A7',([r,]) := A^r?), (1-4) 

is a group morphism where 17 is any continuously dijferentiable representative of [17]. 

The practical implication of this theorem for calculations is as follows: given a model H : 
Q — > A, a value for the Fermi energy E and a deformation 17, all we need to figure out is 
the equivalence class [17] e tiiCQe). Then we are free to use any loop 17' in the equivalence 
class [17] to compute the polarization. Moreover, we get an immediate criterion for the 
triviality of the polarization: 

Corollary 1.2 A necessary condition for AV(,r]) ^ is ti^CQ^) 7^ {0}. 

If TiiCQg) = {0}, then all deformations supported in Q^ are equivalent to the case of no 
deformation, and hence A.V{r]) = 0. It is in this sense that regions in Q where the spectral 
gap at E closes create the non-trivial topology necessary for a non-trivial piezo effect. 

The linearity of AT-"* implies it sends commutators of loops to 0, and thus ( |1.4| l is also 
well-defined as a map from the abelianization of tt^CQ^) to Z''. The abeHanization, how- 
ever, corresponds to the homology group H^(Qe) | jHat02j . This distinction, which is only 
important in the 'exotic' cases of non-commutative ttjCQ^), shows that the piezoelectricity 
depends more properly on the homological properties of Q^. 

1.2 Tight-binding models for piezoelectricity in graphene 

After answering (Ql), let us turn our attention to graphene and the second question. To 



the best of our knowledge, with the exception of the Rice-Mele model in d = 1 [ RM82 



OMN04 1, there are no other concrete models in d > 1 for which the polarization has been 
calculated exactly. Our framework allows for the evaluation of Chern numbers (including 
the polarization) for so-called two-band hamiltonians (cf. Section |4]); in particular, we 
will consider the uniaxial strain model in d = 2. This simple model incorporates all the 
hallmarks of piezoelectrical modifications of graphene proposed by theoretical physicists 
[ ,0R12J , and we show that it allows for deformations which have non- trivial polarizations. 
Let us quickly recount the basics of the crystal structure of graphene: it is an essentially 
two-dimensional material consisting of a single layer of graphite. The carbon atoms are 
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Figure 1: The honeycomb lattice as the superposition of two triangular lattices F with 
atoms of t3^e (white) and type 1 (black) . Every unit cell V contains a pair of 
sites one of t)qDe and the other of type 1. The green arrows 5q, 5^, 5^ connect 
sites of type with the three nearest-neighbor sites of type 1. The blue arrows 
±Yi,±Y2,:i^Y3 connect a given site with the six second-nearest-neighbor sites. 
ip'^^^ (resp. i/)'-^'') denotes the component of the wave function for sites of type 
(resp. 1) in the cell located in j- 



arranged in a honeycomb lattice (cf Figure [T]) which is obtained by the juxtaposition of 
two triangular lattices F c^Z'^ generated by the fundamental vectors 



ri 



f(3, + y3), 



72 



; (3,-73). 



Here a ^ 1.42 A is the distance between two carbon atoms. The vectors 

5o = a(l,0), 5i = |(-l, + V^), 52=|(-l,-^/3), 

connect nearest-neighbor sites belonging to different sublattices. The presence of two 
atoms per unit cell can be described by an internal degree of freedom usually referred to 
as isospin. Hence, if we ignore the electron's spin, the relevant Hilbert space is l^(r)®C^, 
the space of square summable sequences on the lattice F with an internal C^ isospin de- 
gree of freedom. 
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The simplest model which includes only nearest-neighbor hopping is defined by the 
hamiltonian 

V.l«Hr) + qiSi+q2S2 J 

where for j = 1, 2 the operators 

are shifts by jj and the q^ e K are amplitudes which quantify the hopping to nearest neigh- 
bors located in adjacent unit cells. We have fixed the hopping amplitude corresponding 
to shifts by 5q to 1 by fixing a suitable energy scale. The isotropic case qi = q2 = 1 
is the standard tight-binding model for graphene (up to a rescaling in energy of order 
^ — 2.8eV). In configurations (qi,q2) close to (1, 1), the band spectrum of ( |1.5| ) has two 
conical intersections and no spectral gaps. 

In this framework, we assume the net effect of applied strains is captured as a change of 
hopping parameters (qi,q2)- The range of validity of this approximation has been studied 



extensively [ PCP09 RPP+09| , and in our units, it suffices to consider the q^ in the range 
[0,2]. The dependence of the spectrum of T[q-^,q2) on the hopping parameters is well- 
known fHRN-t-OGl (cf. Figure^]): o"(r(qi,q2)) is symmetric around the zero energy, and 
thus the relevant Fermi energy E = lies directly where the spectral gap will open. If we 
restrict ourselves to positive hopping parameters, then the part of the parameter space Qg 
where the gap is open is comprised of three disjoint simply connected components. Thus, 
TTiCQo) = {0} for each component and according to Corollary |1.2[ the piezoelectric effect 
has to be absent. The presence or absence of topological invariants is closely related to 



the symmetries of a system [[AZ97 SRF-l-08 1 : absence of inversion symmetry is a necessary 



condition for a material to be piezoelectric. However, r(qi,q2) has an inherent inversion 
symmetry. Let p be the unitary operator defined by [pip)y := i/'-y Then pSjP = s* holds 
and a simple computation yields that if we tensor p with the Pauli matrix o-j, we obtain 
an inversion symmetry of r(qi,q2), 

[Tiqi,q2) , P®(^i] =0. 

In other words, graphene is not intrinsically piezoelectric. To have any hopes of seeing 
piezoelectric effects, graphene needs to be modified in such a way as to break its inherent 
inversion symmetry. One potential way to achieve this is to adsorb atoms on one side 
of the graphene sheet (e. g. hydrogen, lithium, potassium or fluorine); the piezoelectric 
effect of modified graphene is then expected to be comparable to that of 3d piezoelectric 
materials [|0R12| . The simplest way to capture this breaking of inversion symmetry in the 
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Figure 2: (a) Representation of the parameter space [0,2]^ for the hamiltonian T{(li,q2) 
given in ( |1.5| ) at Fermi energy £ = 0. For values of the parameters in the (closed) 
red region the system is gapless. The gapped part of the parameter space is made 
by the three disjoint triangular regions in white and each of this region is simply 
connected, (b) Representation of the parameter space Q := [0,2]^ x [—1, 1] for 
the hamiltonian H(qi,q2,q3) given in ( |1.6| ) at Fermi energy £ = 0. The extra 
dimension q^ given by the stagger perturbation allows the gapped parameter 
space Qo (i. e. Q minus the red region) to be path-connected. 



model is to add a stagger potential to T, i. e. to consider the uniaxial strain hamiltonian 



W(qi>q2,q3):=7'(qi,q2) + q3 



+1/ 



^(r) 







-1 





ecD 



(1.6) 



instead. We take the parameter space to be Q = [0,2]^ x [—1,-1-1]. Now the gapped pa- 
rameter region Qg is arcwise connected and has a non-trivial fundamental group 7ii(Qo) — 
Z^ (cf Proposition 5.21. Hence, any loop 77 : §^ — > Qq can be continuously deformed 



into a loop which winds n^ times around (1,0,0) and rij times around (0, 1,0) (cf Fig- 
ure |2](b)), and we get 



AVin) = "1 A7'(t}i) + n2 AVini) 



(1.7) 



where 17^ are the loops indicated in Figure[2|(b). In order to prove that this model supports 
non-trivial piezo effects, we need to show AVirjj) 7^ 0. In Sections l4| and [sl we develop a 
technique in the spirit of [Koh85| which allows us to compute the AV[r]j) (and all other 
Chern numbers) explicitly. 
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Theorem 1.3 (Piezoelectric effect in the uniaxial strain model) There are periodic de- 
formations 17 of ( |1.6| ) such that AV{r\) 7^ 0. More precisely, let 



ViM := (1 + e cost,0, — e sint), ViM •= (0, 1 + e cost,— e sint), 

be the two generators o/ttj^CQq) ^ Z^ for some e e (0, 1) andHj^^t) the periodic deformation 
of the graphene Hamiltonian ( |1.6[ ) along r\j. Then AV{r\i) = (1,0) and AV{r\2) = (0, 1), 
and thus AV{r]) = {n^, n2)for [17] = n^ [tJj] + rij [rjj]. 

For details of the calculations, we refer the interested reader to Section [s] 

1.3 Stability under weal< perturbations 

The uniaxial strain model above discussed is based on two simplifications: independence 
of electrons and absence of impurities. A more realistic model should include those aspects 
as well. Mathematically, we can include these effects by adding a potential V to the uni- 
axial strain hamiltonian ( |1.6P , 

H;i(q):=H(q) + AV. (1.8) 

We assume V e ^ is bounded (||V||^ = 1 for simplicity); The perturbation can describe 
interactions between electrons in a mean-field approximation (periodic potential) as well 
as the effect of impurities (Anderson-t3^e potential). The parameter A describes the 
strength of the perturbation. Standard perturbation theory says that if the distance be- 



tween (7(H(q)) and E is greater than g > A, then E ^ o"(H;^(q)) [Kat95|. If Q^ denotes 



the gapped parameter space for the unperturbed hamiltonian H(q) then the set 

Q£_g := {q e Q^ | dist(CT(H(q)),£) > g} (1.9) 

is certainly contained in the gapped parameter space of the perturbed hamiltonian Hp^[q). 
In the weak perturbation regime A e [0, A^] with A^ < g, the space Q^ is a deformation 
retract of the space Q^ g and so the two have same homotopic type [ Hat02 1 . Given a loop 
17 in Qg g one can define two periodic time-dependent and gapped operators H^[t) and 
H}^.^[t) according to the prescription ( |1.2p . The Fermi projections associated to these two 
hamiltonians are homotopic in the algebra A (Proposition |2.4p . As a consequence of the 
homotopic invariance of the King-Smith-Vanderbilt formula ( |1.1| ) [ |ST13 Corollary 2] one 



deduces that H^(^t) and H;^^^(^t) produce the same polarization vector [ ST13 Corollary 3] . 
This fact can be stated as follows: 

Theorem 1.4 Piezoelectric effects persist under weak perturbations. 

This result applies directly to the case of the strained graphene (see Figure [3JI. 
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Figure 3: Representation of the gapped parameter space for the model ( |1.6D perturbed by 
a bounded potential AV in the regime of a weak perturbation A < 1. The topol- 
ogy of this space agrees with the topology of the unperturbed gapped parameter 
space. 



1.4 Organization of the paper 

First, in Section |2J we will reformulate the problem in an algebraic language. Among 
other things, this allows us to include effects of weak disorder just as in [ST13J . Next, in 
Section|3]we sketch the derivation of the King-Smith-Vanderbilt formula dl.ip and discuss 
its topological nature. In particular, we connect ( |l.ip to the topology generated by the 
gap at E in parameter space. Then, two-band systems are discussed in Section|4} we show 
how to exploit the fact that the Bloch bundle can be written as the pull back of a reference 
bundle over the 2m-sphere (the Hopf bundle). Finally, we compute the polarization for 
the uniaxial strain model in Section [S] 
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2 Algebras of observables 

This section serves to introduce the main features of the models in which we are inter- 
ested. The use of a C*-algebraic approach allows us to formulate our results both, for 
periodic and random models simultaneously. 

In the following we will deal only with lattice modeb, i. e. with systems with an under- 
Ijdng geometry described by d -dimensional lattices 

r := {r e K'' I r = X -=1"; Tj, " = ("i, . . . , n^) e Z'*} 

generated by d linearly independent basis vectors {Yi,..-,Yd} ^^d set \y\ := 2,=iI"j(t)I- 
We will consistently use the notation n^j) •= ("iCt), •••,nd(T)) ^ ^'^ for the vector of 
coefficients which express 7 e F in terms of the basis {ji, . . . , 7^}. 

In order to include internal discrete degrees of freedom like spin and isospin we will 
usually consider tensorized objects of the form 

2l:=*B ® Mat,(C) (2.1) 

where ® is any complex normed (or locally convex) algebra and Mat^CC) is the algebra of 
the r X r matrices with complex entries. The fact that Mat^CC) is a finite dimensional al- 
gebra (hence nuclear) implies that the topological tensor product ( 12.1| ) is uniquely defined 
without ambiguities [ |Tre67[ . Moreover the following identification 

*» ® Mat,(C) ^ Mat,(23) 

will be tacitly used when convenient. 

2.1 The algebra of periodic operators 

There is a canonical way to have a unitary representation of the lattice F in terms of shift 
operators on £^(F): to each generating vector y^, j = 1, . . . , d, we define 
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2 Algebras of observables 



The shifts commute amongst each other and we can use multi-index notation to define 
the group action 

d 

B-.r ^Bil\n), r^s,:=Y\5f'\ (2.2) 

;=i 

Starting from the algebra of finite linear combinations of shifts 

6fi„:={ae6(£2(r))|3JVeNo : a = S.^i^^a^s,}, 

we define the shift algebra & as the completion of &f^^ with respect to the operator norm 
on £^(r). Including also the internal degrees of freedom one defines on the spinorial 
Hilbert space 

-H — f^Cr)®^ 

the tensorized algebra of periodic operators 

%^,:=e®MatXC)^MatX&). (2.3) 

To avoid confusion between elements of the 'abstract' Brillouin algebra Ap^^ (to be defined 
in Section 2.2 1 and its representation 2lpgi. we denote elements of Slpej. with a hat. The C*- 
algebra 2lpgf admits a differential structure and an integration. The position observable 
X := (xi,...,Xjj) is the vector-valued (unbounded) operator defined component-wise on 
i\n by 



(xja/))(r) = n/r)^(r)- (2.4) 

The position operators let us define (spatial) derivations on 2lper by 

VjA:=i[A,Xj®l,], j = l,...,d. (2.5) 

Clearly, the V^'s are unbounded operators as can be readily seen from 

Vj(s^®M) = -in^(r)s^®M (2.6) 

for every M e Mat^CC). Nevertheless, the Vj's are initially defined on the dense subal- 
gebra 6^^ ® Mat^CC) and then extended to their natural domain. As can be checked by 
explicit computation, these derivations are symmetric, (V^A)* = Vj(A*), and satisfy the 
Leibnitz rule Vj^AB) = Vj(A)B +A'V0). It is also easy to check that V^ and V^ com- 
mute. The gradient operator V := (Vj, . . . , Vj) is closable on a common natural domain. 
Let us define the norms 

ll^llp ■■= S W^^AWm ^2.7) 

laNp 
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2. 1 The algebra of periodic operators 



where for any a = (a^ . . . a^j) e Ng we used a multi index notation V" := V"' • • • V^'' with 
|a| = |ai| + ... + |a^|. The completion of the dense subalgebra 6(5^® Mat^ (C) with respect 
to IHIp leads to the Banach-* algebra oi p-times dijferentiable operators C^(2lper) c 2lper- 
From the definition it follows that C+HSlper) '^ C^CSlpg^) for allp e Nq with the convention 
that C°(2tper) = Slper- The algebra C^(2tper) is the natural domain for V while the Frechet- 
* algebra C°°(2lper) := nneN^^^^^per) (endowed with the inductive limit topology) is an 
invariant domain of V. 

Remark 2.1 All the algebras CCSlpgJ are stable under holomorphic and continuous func- 
tional calculus. Moreover, one can prove that if H = H* e CCSlper) and / e C+H*) 
hold, then /(H) e 21 r, defined through continuous functional calculus, is in fact an el- 



ement of C^(2lper) [ BEJ84| Lemma 3.2]. This does not only hold for Ar,^^, but also all 



non-commutative algebras which we will work with in this paper. 

The second relevant structure, the integration, is defined on 2lper by the so-called trace- 
per-unit-volume T. Let 5^ ® e^ be the canonical basis of H where {ej, . . . , e,.} denotes the 
canonical basis of C and 5^ := (5^^')^'gp is the £^(r) normalized sequence with only one 
non-zero entry at the label 7. Then the trace-per-unit-volume of A e 2lper is given by 

r 

nA) := Y,{5o®epA5o®e^)^. (2.8) 



This map is a *-linear functional T : 2lpgr — > C with the trace property [ ST13 Lemma 1] . 



The name trace-per-unit-volume is justified by the following observation: Let {Fnlngfsj be 



any F0lner sequence [F0155 1 of bounded subsets of the lattice F such that F„ c F„+i and 
F„ /■ F and denote with |F„| the cardinality of the finite set F^. If one introduces the 
orthogonal projections 

rer„ 
and if one observes that Xy = (.^y ® 1) Xo (^^y ® Ir)*- one can verify from ( |2.8| l that 

r(A) = Tr„(fo^Io) = T|7-|Tr«(zr„^fr,.) =^lm r^Jr^ix^^Ax^J- 

The calculation uses in a crucial way the cyclicity of the trace and the translation invari- 
ance of A e 2lper, namely [s^ ® 1^ , A] = for all 7 e F. The last term in the above equality 



provides the usual definition of the trace-per-unit-volume (see [ Ves08 1 for a general re 



view) and the equality does not depend on the particular choice of F0lner sequence. 



13 



2 Algebras of observables 



2.2 The Brillouin algebra for periodic lattice systems 



According to ( |2.3| ) the non-commutative part of the algebra Slp^^ comes entirely from the 
matricial factor Mat^CC) since the algebra 6 is commutative. This last observation allows 



us to use the Gelfand-Naimark theorem [H6r90|: it estabhshes a C* -algebra isomorphism 



between 6 and the C*-algebra of continuous functions C(Spec(6)) where the algebraic 
spectrum Spec(6) is a compact topological Hausdorff space. Since the C*-algebra © is 



generated by the Sj, we can use [DP 12 Proposition 5.5] to write 



Spec(6) = Y\ (^^Sj) = S^ X • • ■ X §1 = T'*. 

j=i 

The Gelfand isomorphism ig : & — > CCT'') is uniquely defined by its action on the genera- 
tors Sj -^ e"'''', and it maps Xiver '-r^r '^^'■^ ^^^ corresponding trigonometric polynomial. 
Reversing the direction of the isomorphism yields a representation of the C* -algebra 
C(T'') onto the C*-algebra of operators 6 on the Hilbert space £^(r). This representation 
extends after tensoring with the matricial part Mat^CC). More precisely, let us define the 
periodic Brillouin algebra 

Ap^r ■■= CCT'') ® Mat,(C). (2.9) 

Using the identifications A^ei - MaX^{C(J'^)) and Slper ^ Matr(6), and the Gelfand iso- 
morphism ig^ for each component, we define a faithful representation of the periodic 
Brillouin algebra on the algebra of periodic operators, 

7Zp,,:Ap,,^%,,^B(ni (2.10) 

As we shall explain in the next subsection, this point of view extends in a natural way to 
the case of random operators. The representation iZp^^ can be concretely realized through 
the Fourier transform 

which is a unitary map J" : £^(r) -^ L^^T^) between Hilbert spaces. A simple computation 
shows that J^Sj F* = e"'''^ where the right-hand side must be interpreted as a multipli- 
cation operator on L^(T''). This means that the Gelfand isomorphism i^^ is unitarily 
implemented by i^^^f) = J^ f J- ^& ior each continuous function / in C(T''). Tensoriz- 
ing the Fourier transform by the identity matrix 1^ one obtains also a unitary description 
of the representation ( |2.10| ). More precisely for each continuous matrix-valued function 
A e ylpgf one verifies that 

np^,{A)^{F®\,rA{F®\,)=:A&%^,. 
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2.2 The Brillouin algebra for periodic lattice systems 



The fact that the algebra of operators Slpg^ is just a faithful representation of the algebra 
Apf.^ allows us to to investigate spectral and dynamical aspects directly in the algebra 
Apf.^. The first advantage concerns the calculation of the spectrum. The above relation also 
means that A e Ap^^ andA= Tip^^{A) e Stpej. are unitarily equivalent, and thus cr(A) = cr{A). 
Now the spectrum of ^4 is easily accessible since it acts as a matrix- valued multiplication 
operator on the fibered Hilbert space L^^T^) ® C. Of particular interest is the case of a 
selfadjoint H = H* & Slp^^ . In this case for each k&T^ the operator H(fc) is a symmetric 
r X r matrix with eigenvalues £i(fc) ^ ... ^ E,.[k)- The functions k -^ Ej[k) are called 



energy bands. As a standard result [RS78| we have a complete characterization of the 
spectrum: 

Lemma 2.2 The spectrum of any selfadjoint H =H* & 2tper has empty singularly continuous 
components and consists of closed intervals, 

o-(H) = /i U...U/r 

where Ij := U^.t^^;^^)} = [Ef\Ef''^]. 

Also the differential structure and the integration defined on the operators algebra Slp^^ 
have a counterpart on the level of the algebra Ap^^- First of all, since any A in ^p^j. can be 
seen as a map from the manifold T'' to the normed algebra Mat^CC) the d partial deriva- 
tives d^A are well defined (assuming A is sufficiently regular). A formal computation on 
linear combinations of generators and elementary tensor products leads to the relation 

VrC V) = i [^perW ,Xj®l,]= V^ (7Tp,,(A)) . (2.11) 

With the identification of the notation d^, = V^ we can rewrite the equation ( |2.11D as 

7Tper°V; = V,0 7rp„ (2.12) 

which means that the representation Tip^j. intertwines with the gradient V. The above 
relation is well defined on a maximal domain. Let us introduce the regular subalgebras of 

A ■ 

CPIiAper) = CP(T^)®Mat,[C) ^ Mat,{CP(T^)), p = 0, l,...,oo. 
Then it is straightforward to check that Tipg^ defines faithful *-algebra map 

This together with equation ( |2.14| ) establishes that C^^Ap^^) ^ -^per is the natural domain 
for the gradient V and C°°(^per) c Ap^^ is an invariant domain. 
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2 Algebras of observables 



On Aper — Mat^ [C(l"^)) the integration T involves a bona fide integral and the trace, 

'^^^)-= I 7:r^Tre(A(fc)). (2.13) 

Jjd [271 j" 

Here dfc is normaHzed such that f ^ dfc = [2tzY. Also in this case one can verify (first on 
a dense subalgebra) the intertwining relation 

T = ToTip^, (2.14) 

between ( |2.13| ) and the trace-per-unit- volume ( |2.8| ) . In particular all properties listed in 



[ |ST13| Lemma 1] hold true also for ( |2.13D . 



Remark 2.3 We stress that the use of the symbols V^ and T both for the algebra Ap^v 
and its realization 21per has the great advantage of allowing us to write formulas like 
( |1.1D independently of the specific realization in a given algebra. Moreover, the risk of 
confusion caused by this abuse of notation is minimal and, when necessary, the reference 
to the algebra will be mentioned explicitly. 

2.3 Covariant families of random operators 

A random system on the lattice F is described by hamiltonians of the form 

where H is a periodic operator, i. e. an element of 2lpg„ V^ is a bounded operator on 
v. which depends on a random parameter co and a coupling constant A > 0. A typical 
example is the Anderson potential 

V„ := 2^,15^) (5,1 ®1, 

with CO, e [0, 1]. The collection co := (w,),gr defines a configuration of the disorder The 
configurations take values on the space n := [0, 1]^ which turns out to be compact if 
endowed with the Tychonofi^ topology. If all the one site configurations w, are distributed 
on the interval [0, 1] according to the same probability measure d/i, one can endow also n 
with the product probability measure dP := x,gp dfx. The group F acts on the topological 
space Q by translations via (t^o)),/ := co,'_,. The measure dP is invariant and ergodic 
with respect to the group action t. Moreover, using the invariance of H with respect to 
the translations Sy® 1^, one can verify the covariance property 

(s,®l,)H„(s,®l,r=H.„, VreF. (2.15) 
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2.4 The non-commutative Brillouin algebra for random lattice systems 

The main features of the Anderson model can be used in order to provide a general def- 
inition of random lattice systems. Following the above example, the randomness can be 
described by a triple (r2,dP, t) where fi is a compact Hausdorff space, dP is a borehan 
probability measure and t is an action of F on fi by homomorphisms. The measure dP 
is required to be invariant end ergodic with respect to t. Associated with this structure 
we can consider family of bounded operators (Aco)oj€n ^ BCH) such that: (i) the map 
CO -^ A ^ is strongly continuous and (ii) the covariance property ( |2.15P holds true. We 
refer to such a (A„)„gfj as a covariant family of random operators. We stress that in view 
of ( |2.15D periodic operators can be identified with constant covariant operators, namely 
with a random family such that j4„ = A for (almost) all co G fi. 

Instead of a particular realization, one studies the covariant family of random operators. 
Many of their important properties are in fact deterministic, e. g. spectrum and spectral 



components [ PasSO KS80| 



a(A^) = S, P-a. e. wen. 



and the P-a. s. existence of the trace-per-unit-volume [Bel86| 



nAj:= hm -^Tr^(xr,.^Ir„)= dVTr^(xoA^ Xo)- (2-16) 

Both of the above properties are consequences of the Birkhoff ergodic theorem. 

Also the notion of derivative extends to random families of operators as a P-a. s. prop- 
erty and gives rise to p -times differentiable and smooth covariant families of random 
operators. 

2.4 The non-commutative Brillouin algebra for random lattice systems 

The aim of this section is to construct an 'abstract' C*-algebra which encodes all the topo- 
logical and geometrical characteristics of the set of covariant random operators. This 
algebra can be thought of as a generalization of the C*-algebra of periodic observables 



Apgj described in Section 2.2[ This construction has been developed by Bellissard during 
the 1980's [ |Bel86[[Bel88 |. 



Given a covariant family (A^)^^q, the main idea is to consider ^^(w, y) := {5q,A^ ^r)« 
as r X r matrix-valued symbols for covariant operator families and to construct a C* -algebra 
out of these symbols, given by an adequate crossed product. First one endows the topo- 
logical vector space C^. (fi x F, Mat^CC)) of continuous functions with compact support on 
n X F and values in Mat^CC) with a *-algebra structure: 

AB[co,r) := 2A(w,/)B(T_,,w,r-r')> (2.17) 

r'er 

A%co,r) := A(T_,w,-r)*. (2.18) 
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2 Algebras of observables 



For any w e fi, a representation of this *-algebra on V. is given by 

(7r„(A)*)(r) := Y,^iT_^co,r' - r) *(/) (2.19) 

r'er 

where * e H and ^(y-) e C for all 7 e T. From ( |2.19D it follows that the different 
representations 7r„ are related by the covariance relation 

(5^ ® 1,) Tl^iA) (S^ ® 1,)* = 7I^^„(A) 

and are strongly continuous in o). With the notation ti^(A) =: A^ one can see that the 
family of representations n^ sends elements of the *-algebra C^[n x r.Mat^CC)) to co- 
variant families of random operators. 

If we complete C^{n x r,Matr(C)) with respect to the C*-norm 

\\A\\ :=sup||7r„(A)||g^^j, 

we obtain the (non-commutative) Brillouin algebra 

A := (C(n) x: r) Mat,(C) ^ Mat, (C(n) x: r) (2.20) 

where C(n) x: F denotes the cossed-product C*-algebra [ |Wil07[ . Note that since n is com- 
pact, the algebra C(r2) is unital. Hence, the crossed product 

Ao := (C x: F) ® Mat,(C) =^ Mat, (C x: F) (2.21) 

is a C*-subalgebra of A which consists of the elements that are independent of co. 

The algebra A carries a differential structure and an integration. The first is given by a 
gradient V := (Vj, . . . , V^ defined by 

{VjA){co,r) := injiY)Aico,r). 

The domain of this gradient is the subalgebra C^[A) which is the completion of the dense 
algebra C(.(f2 x F,Mat,(C)) with respect to the Banach norm ||-||i given by a formula 
analogous to ( |2.7| ). More generally, with the usual procedure, one can define also the 
algebras of p -times differentiable elements C^^A) and the algebra of smooth elements 
C°°[A) which is an invariant domain of V. A simple computation provides 

n^iVjA) = i [n^(,A),x^ ® 1,] = V^{7ip,,(,A)). (2.22) 

which means that the representations n^ intertwine with the gradient V, 

7T„°Vj = Vj0 7T„. (2.23) 
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2.5 Perturbations of periodic operators 



An integration is defined on A by the tracial states 

r(A):= dPTre(A(w,0)). (2.24) 



A comparison between the definitions ( |2.19P and ( |2.16[ ) provides also in this case the 
P-a. s. intertwining relation 

r=r°7i„. (2.25) 

2.5 Perturbations of periodic operators 

To consider perturbed periodic operators, we start by showing how to identify Aq (z A 
made up of elements independent of co (cf. equation ( |2.21| )) with Ap^^: Equation ( |2.19| ) 
says that every n^ maps elements of Aq to bounded operators in V. which commute with 
translations Sj ® 1^, i. e. to periodic operators. Periodic operators in turn are represented 
faithfully by 71^^^. Following this reasoning, we see that 

establishes an isomorphism between Ap^.^ and Aq, and that this isomorphism does not 
depend on co. 

Starting with an hamiltonian H = H* & Ap^^ ^ -^ which has a spectral gap at E (gap 
condition), 

dist(£,a(H)) ^ g>0, 

we can perturb it by a bounded periodic potential in such a way that Hp^ = H + A Vp^f is still 
an element of ^p^^ or by a bounded random covariant potential so that Hx= H + AV„ is 
now an element of ^. In both cases H^ converges to the unperturbed periodic hamiltonian 
H in ^ as A ^ 0. Consequently, standard perturbation theory in the sense of Kato [ Kat95| 
guarantees the persistence of the gap as long as the perturbation is not too strong: there 
exists A^ < f such that 

* 2 

dist(£,CT(Hj)^| 

holds for all A e [0, A^] . This means, we can define the Fermi projection P;^ := 1(-oo,£](^a) 
for all A e [0, A^]. By standard results, P^ is also an element of A (cf Remark 2.11; 
furthermore, if the perturbation is periodic, then the Fermi projection is also in Ap^.^- 
Then the continuity of A -^ P;i can also be interpreted in the following way: 

Proposition 2.4 Under the conditions listed above, [0, A^] 9A'-^P;^eylisa homotopy. 
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3 The King-Smith-Vanderbilt formula 



Proof The continuity of the family of bounded operators Hp^ in A and the resolvent iden- 
tity imply the local continuity of the resolvents [Hx — z)~^. The gap allows us to write 
the Fermi projection as a Cauchy integral involving a contour that can be chosen indepen- 
dently of A. Hence, A -^ P;^ is continuous or, in other words, it is a homotopy □ 

The importance of this result resides in the fact that all physical quantities which depend 
only on the homotopy class of the spectral projection can be computed in the limit of zero 
disorder. This meta result is generally known as stability under weak perturbation. 

3 The King-Smith-Vanderbilt formula 

The model hamiltonians H : Q — > A we are interested in are parametrized by a parameter 
space Q. The latter is always a suitable path-connected C^-submanifold of K'^ where 
notions such as taking derivatives are tacitly inherited from R'^]^ yl is a *-algebra. We 
shall always make the following technical 

Assumption 3.1 We assume H -.Q czR^ — > C^(.A) is a C map with p ^ 3 taking values 
in the self adjoint elements ofC^(A). 

Let us pick a Fermi energy E and consider the gapped parameter set Qg defined as in 
the introduction. Then the relation < \1.2\ i mediates between loops 17 e C(S^,Qe) and 
time-dependent, T-periodic hamiltonians H^[t). Note that if Qg has several connected 
components, then 71 ^(Q^) is the direct sum of the fundamental groups of the connected 
components. 

3.1 Derivation from first principles 

To give a self-contained presentation, we will sketch the derivation of ( |l.ip . The assump- 
tion H^[t) e C^iA) implies that the current operator 

;^(t):=VH^(t) 

is bounded. 

The dynamical polarization is the expectation value of the charge transported over one 
period, and a quick computation | jST13[ Proposition 4] yields 

APdyn(i7):= dtr(P(OVH^(t)) 
Jo 

• r 

dtr(p(t)[5,p(t),v/(t)]). 



/o 



^Indeed, one may also consider parameter spaces which are Riemannian manifolds, but for our intents and 
purposes, this is not necessary. 
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3.2 The polarization as a topological quantity 



where P[t) is the solution to the Liouville equation with initial state P^(0) as given by 
equation ( |1.3| ). 

Assuming the deformation is sufficiently slow and regular, we can approximate AT'jjynCrj) 
with AV(,ri) by replacing the time-evolved projection P[t) with the Fermi projection P^^^t): 
if t -^ H^[t) is C as a T-periodic map from K to C^[A), then LSTIS. Theorem 1] states 
that the error is of (p — 2)th order in the adiabatic parameter e, 

3.2 The polarization as a topological quantity 



A second, and for our purposes equally important result, [ST13 Corollary 2], says that 
AV is invariant under C^-homotopies oi projections which gives us leeway in how to calcu- 
late AV. The remainder of Section[3]serves to show that instead of looking at homotopies 
oi projections, it suffices to look at homotopies in parameter space. 

Lemma 3.2 If 17 and 17' are in the same equivalence class of the p-regular homotopy group 
TTjCQg), then also the Fermi projections P^ and P^i are C^ -homotopic. 

Proof Let 17,17' e CP(S^,Qe) be two loops with [17]^ = [17']^ e tIjCQ^). By definition, 
there exists a C^ -homotopy 

which connects 17 = A(0) with 17' = A(l). The continuity of (s, t) -^ Hj^r^^^t) ensures the 



inner and outer continuity of a [H jy(^^-j[t)) (see e. g. [AMPIO Corollary 2.6]); moreover, 



the resolvent (s, t) -^ [Hj^^^^^^t) — z) inherits the C regularity of q -^ H(q) and A(s). 
Hence, writing P^j-j^Ct) as a Cauchy integral, we see that the map (s, t) -^ P^jj-iCt) is also 
CP. 
Now we cover [0, 1] with finitely many open intervals {V^} such that 

sup||PA(,)(t)-PA(,,)(t)||<l (3.1) 

holds for all s,s' e V„. 

Let us initially assume ( |3.lD holds for all s,s' e [0, 1]. Then this condition implies the 
existence of a family of unitaries [ |Kat95| equation (4.38)] 

[/(s; t) := (^P^^,^^t)P^^t) + (1 - PA(,3(t)) (1 - P^(,t))) (1 - (PA(,)(t) - P,(t))')"'^' 

(3.2) 

which intertwines P^^iit) and 

Pa(s)(0= [/(s;t)P^(t)[/(s;t)*. (3.3) 
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3 The King-Smith-Vanderbilt formula 



This unitary [/(s; t) inherits the C regularity from Pa(s)(0 ^nd P^(t). 

Let us return to the general case: If we need several open intervals {V„} to cover [0, 1] 
so that ( 13.1| ) is satisfied on each of them, the above procedure yields a collection of ho- 
motopies {Pa^s; t)} where each of the P„(-; t) is only defined on V„. Gluing these homo- 
topies together yields a homotopy P[s; t) connecting P^[t) and P^^'^t) which is C^ almost 
everywhere, but only continuous at the gluing points. However, we can invoke [|BT82| 
Corollary 17.8.1] and modify the homotopy P(-; t) to make it C^ everywhere. □ 



According to [ ST13 Section 3], the quantity AVirj) given by ( |1.1D (for a fixed 17) is a two- 
cocycle on the extended C*-algebra A := C(§^) ® A endowed with the extended gradient 
V := (i3f, V) and the extended trace T := J dtT(0- Such an object provides the non- 
commutative analog of the Chern invariant in the spirit of non-commutative dif ferential 



calculu s. More precisely, T can be seen a map between the i<^-group Kq(^A) and Z [ Con94 



VFGOl]: T applied to a projection in A yields an integer, and if two projections are Kq- 



equivalent, they are mapped to the same integer. On the other hand. Lemma [3^ says that 
for 17 and 17' from the same equivalence class of the p-regular homotopy group tIjCQ^), 
the Fermi projections P^ and P^, are homotopic in A, and so in the same JCo"Class. This 
leads to the following result: 

Proposition 3.3 For any integer p ^ 1 equation ( |1.1| ) induces a group morphism 

AV::CP[S\Qe)^Z'', [rj]p-A7'f([r,]p) := AV[r^l 

The linearity of AV^ follows directly from the linearity of AV; also APf ([0]^) = is 
immediate, because B^Pq = for Fermi projection associated to the constant loop 0. 

3.3 The periodic case: thie Blochi bundle 

In the periodic case, the polarization AV can be seen as Chern numbers associated to the 
so-called Block bundle, and the C^-diffeotopy invariance is but one consequence of this 



fact. Using the notation of Section 2.2 the Fermi projection P^^^t) e Ap^r ^ Mat^ (C(T'')) 
can be viewed as a family of projections PJk, t) on C indexed by (fc, t) e T''. Then in 



analogy to [ Nen83 ; Pan07 DLllJ , the disjoint union 



£bW-= U ranP^(fc,t) (3.4) 

defines the so-called Block bundle. The continuity of (fc, t) -^ Pj,(fc, t) directly implies 

7r:£B(r,)^T''+i 
is a vector bundle. This construction is nothing but a particular application of the Serre- 



Swan tkeorem [Swa62|. Since homotopic loops define i<^o-equivalent projections in the 
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3.4 The topology of the parameter space 



module C(§i) ® Matr(C(T'')) c^ Mat,(C(T''+i)) (cf. Lemma 3.21 and isomorphic vector 
bundles arise from i<^o -equivalent projections [ |Hus66 |, an immediate consequence is 

Proposition 3.4 The Block bundle n : ^bC^) — * T''"'"^ depends only on the homotopy class 
[17] e ttiCQeX i. e. if 17 ~ 17', then the bundles £^[ri) and E^irj') are isomorphic. In particu- 
lar, all Chern classes agree, c„(£'b(i7)) = c„(£'b(i7'))- 



We will nov\? connect AV to the Chern class: the homology group HjCT "*" ) = Z 



■id(d+l) 



[ DLll Section VC] is generated by the 2-tori 

i2 



T^ 



J.n 



T 



n,d + l 



{(fc,OeT'*+i I {k„...,l^,...,])^,...,k„t)=* where *eT''-i}, 
{(fc,OeT'*+i I ik„...,)^,...,k„/)=* where *eT''-i}, 



where * e T'*~^ can be chosen arbitrarily for each ^ j < n ^ d separately. Here, 
(fcj, . . . ,^, . . . , fcd) e T'*~^ means we omit fc„. Fixing the point * corresponds to a choice 
of embedding T^^ ^ T''+\ and how we identify T^^ with [-Ti,+7zf ^ [-7r,+n;)''+\ 
Clearly, this choice has no bearing on the resulting first Chern numbers 



(3.5) 



Cj,„([i7]):= Ci(5B(r?))eZ, 



and in calculations some choices are more convenient than others. If we set C^ „([i7]) := 
— C„j([i7]) for n < j, define Cjj([tj]) := and arrange 



r 



C([rj]):= 



ObCLt?]) 



V -AVMr)]) 



-AVMr]]) 



+APi([rj]) A 







in an antisymmetric matrix, then the adiabatic polarization AV^liq]) make up the non- 
trivial components of the last column. The remainder f2B([i7]) is comprised of the Chern 
numbers which are relevant for the quantum Hall effect. 

3.4 The topology of the parameter space 



Let us return to the general case including disorder In view of Proposition 3.4 and the 
i<^o -equivalence of homotopic projections, the differentiability condition in Proposition 3.3 
seems superfluous. Indeed, for a wide variety of cases, making this distinction is not 
necessary because ttJCQ^) is isomorphic to tiiCQ^). This is the case, for instance, if Q^ is 
an open subset of R^ , because then, the inclusion map 

i:C°°(,S\Qe')^C(,S\Qe) 
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3 The King-Smith-Vanderbilt formula 



induces an isomorphism i^ : ti^CQ^) — > tijCQ^) between smooth and continuous first 



homotopy groups gBT82[ Corollary 17.8.1]. 
In general, the inclusion map i : C^(S^, Q^) — > C(S^,Qe) induces only a homomorphism 

Very often, it is not necessary to check whether tiJCQ^) :^ tt^CQ^) holds for a variety of 
choice of £, but it suffices to verify n;J(Q) ^ 7ii(Q) instead. 

Proposition 3.5 7iJ(Q) ^ tiiCQ) =^ tiJCQe) ^ ttiCQe) V£ e R 

To prove this, we need a lemma which is interesting in its own right: 

Lemma 3.6 Q^ C.Qis open with respect to the relative topology of Q. 

Proof If Qg = 0, there is nothing to prove. So let qg e Q^ 7^ be arbitrary. Then the 
closedness of o"(H(qo)) as well as the existence of a gap at E implies that for e > small 
enough, the compact interval \_E — <^/2,E + e/2] is fully contained in the gap, i. e. 

C7(H(qo))n[£-^/2,£ + ^/2] =0, a(H(qo))n(-oo,£ + e/2)#0. 

Fundamentally, the continuity of Q 3 q -^ H(q) e A implies the inner and outer continuity 



of the spectra c7(H(q)) (see e. g. [AMPIO, Corollary 2.6]). Then the outer continuity of 
the spectrum a{H[q)) ensures the existence of an open neighborhood Uq^ c Q such that 

a (H[q)) n[E- ^2, E + e/2] = 0. 

holds for all q e U^^ . Moreover, the inner continuity of the spectra guarantees that spec- 
trum does not suddenly collapse, i. e. 

a(H(qo))n(-oo,£ + V2)7^0 

holds on a possibly smaller open neighborhood Uq^. As that also implies E ^ a[H[q)) for 
all q e Uq^, we have in fact Uq^ c Q^. Because qo e Q^ was arbitrary, this shows Q^ is an 
open subset of Q. □ 



Proof (Proposition 3.5 ) Pick any £ e K. Without loss of generality, we may assume 
Qg 7^ 0. As Qg is an open subset of the C^-manifold Q, i^ : ttJCQ^) — > ttjCQ^) is injective. 
It remains to show that i^ is also surjective. 

Let 17 e C(S^,Qg) be a loop. Then we may also interpret 17 as a loop in Q, and because 
Ti\[Q) — T^iCQ), there exist many C^ modifications 17^ to 17 as loops in Q. 

To see that there exist also modifications which lie entirely in Qg, we note that Q^ is 



open by Lemma 3.6 Thus we may pick an open tubular neighborhood U of the graph of 
17 which lies entirely in Qg. Since U can also be seen as an open subset of the C^ -manifold 
Q and 7iJ(Q) ^ tijCQ), we can find a C^ modification 17^ of 17 which lies in K. Then by 
definition [17^] = [97] e tiiCQ^). This concludes the proof. □ 
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Hence, the assumption on Q^ in the fohowing theorem is satisfied in many practical cases: 

Theorem 3.7 (Homotopy-invariance of AV) Assume i^ : 7rJ(Q) — > n:i(Q) is an isomor- 
phism. Then for any £ e K and all continuously dijferentiable loops r\ e C^{S^,A), the 
polarization depends only on the equivalence class [rj\ e tt^^CQ^), i. e. AV induces a group 
homomorphism AT'* : ttjCQ^) — > U^ given hy 

AV,{[r]'Y) = AV{r]) 

where r] is a C^ -representative of [17] e ttjCQ^). 



Proof First of all, Lemma 3.5 states that also n\[Qg) and ti^CQ^) are isomorphic, and we 
will denote the isomorphism also with i^,. The n the composition of the group morphisms 
^ and AVl : n\[Q^) — > Z'' (Proposition 3.3 1 yields yet another group morphism. 



AV.-AVloi^^-.n.iQE^^Z". 



4 Two-band systems 

We have seen that for periodic lattice systems the adiabatic polarization AV(,[r]]) can 
be seen as Chern numbers associated to the Bloch bundle. The goal of this section is to 
provide a technique to compute these Chern numbers with particular interest in models 
where r = 2. 



4.1 Periodic two-band systems 

Aperiodic two-band system is a continuous map H : Q — > ylp„ of the form 

2m+l 

H(/c;q) = ho(fc;q)l2».+ Xi H^'^^^j W.l) 

j=i 

where the functions hj(-;q) e C(T'') are real valued for all q e Q. The 2™ x 2"" matrices 
{T,i,..., S2m+i} sre selfadjoint and provide a (non-degenerate) irreducible representation 
of the complex Clifford algebra Clc(2m)r| Moreover the matrices S^ can be explicitly 
constructed as tensor products of m Pauli matrices. With the above assumptions H(q) := 
n;pgr(H(q)) is a selfadjoint operator on V. for all q e Q. 



^Note that the same set of matrices provides also a degenerate representation for the CUfford algebra Cl(;(2m+1) 



since Si ■ ■ ■ E2„ = (i)'"S2„+i |Lee48 |. 
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4 Two-band systems 



H{k; q) can be diagonalized explicitly: if 

f2m+l \ '1^ 

N(fc;q):= Xi^(fc;q)M (4.2) 

is strictly positive for all fc e T'*, then for this q&Q_ the local gap condition is verified, the 
two eigenprojections 

exist and they satisfy Pj^[k;q)P^[k; q) = 0, P+(fc; q)+P_(fc; q) = Ij". • In view of (H(fc; q)- 
ho(fc; q) l2">) = I^Kfc; q)^ l2"') the matrix H(fc; q) can have at most two distinct eigenval- 
ues. In fact, the orthogonal projections split 

H(fc;q) = H+(fc;q)eH_(fc;q) (4.4) 

into its two spectral components 

H±(fc;q):=P±(/c;q)H(fc;q)P±(/c;q)=(ho±|h|)(fc;q)P±(/c;q). 

For any (fc; q) the rank of Pj-(fc; q) is exactly 2'""^ . This can be checked by simply taking the 
trace of ( |4.3| ) and exploiting that the generators of the Clifford algebra T,j are traceless. 

Remark 4.1 Clearly, for m = 1 all models are two-band systems. However, many other 



models such as the Kane-Mele hamiltonian [KM05 ] are not of the form ( |4.11 l, because they 
also include products of sigma matrices T.j'' := S^^ . . . S^ (with 1 ^ j^ < . . . < j^ ^ 2m 
and p = 2, . . . , 2m — 1). However, if they can be written as 

2m- 1 

H = Ho + 2 XP-'H^, H, := ^ h,^,...,,-^ sjf^'..,^.^, 

p=2 ls:ji<...<jpS;2m 

where Hg is a two-band hamiltonian, then H can be viewed as a perturbation of ( [4.1| ). 
Thus, topological quantities associated to H coincide with those of Hq provided A is small 



enough (cf Section 2.5 1 



4.2 The Bloch bundle of a time-dependent two-band system 

For any loop rj e C(§>^,Qe) we can consider the time-dependent two-band system associ- 
ated to ( |4.1| ) via the prescription ( |1.2| ), i. e. 

2m+l 
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4.2 The Bloch bundle of a time-dependent two-band system 

with h^j{k, t) := h^ (fc; ■)7(2'it/r)) . If for all q e rj the local gap condition |h|(- ; q) > is 
verified (e. g. if there exists an energy E such that 17(8^) c Q^)^ we can define from ( 14.3| ) 
the Fermi projection 

The above formula describes a continuous family of orthogonal projections R/fc, t) on 
C^" indexed by (fc, t) e T'^+^j this family defines the Bloch bundle (cf Section [Jis]) 



£bW-= U ranP^(/c,t) (4.6) 

of rank 2™"^. An efficient way for studying the vector bundle S^^r]) is to introduce a 
spherical parametrization. In fact, an inspection of equation ( |4.5[ ) implies the projection 
depends only on the unit vector u(fc, t) e §^'" with components Uj(fc, t) := 'i.,j('f.O/|h|,|()c,t) 
which suggests to write Pj,(fc, t) in terms of spherical coordinates 

^i(0i,---,02m-i,'/'):=sin0i ... sin02m-2 sin02m-i cos;/? 
y2(0i,...,02m-i.¥'):=sin0i ... sin02m-2 sin02m-i sin(/7 



y/01,...,02m-l,'/'):=Sin0i ... Sin02m+l-j COS02m+2-; (4.7) 

y2m(01.---,02m-l,'/'):=Sin01 COS 02 
^2m+l(01. • • • > ^2m-l. V) " COS ^j 

where the angular coordinates 6i,..., 02m-i range over [0, n] and Lp ranges over [0, 2ti] . 
With the spherical coordinates ( |4.7P one can construct a family of orthogonal projections 
on C^"" given by 

-■ ^ 2m+l A 

Ps2,„(u) := - l2.» - 2 y/u)Ej (4.8) 



with u := [9i,..., 62^-1,^)- The projections ( |4.8| ) define a reference vector bundle n : 
f g2m ^ S^™ over the sphere S^™ with total space 



£g2m := ranPg2,„(u). 
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4 Two-band systems 



This vector bundle is also known as the Hopf bundle in the literature; It is non-trivial with 
Chern class c^(£s2m) ~ {—iy^~^ ^"^~^'' Av where ^2^ is the volume of the sphere S^*" and 



dv e H^^CS^"!) is the volume form (see [Kar87 Exemple 1.27] and [Hat09 Corollary 
4.4]). 
With the help of the inverse functions 



B^jija, t) := arctan 



ip^[k,t) := 2 arctan 



"?j,2m+2-;"Cfcj 



(4.9) 



one obtains a continuous map 

which allows us to write P^ = Pg2m o $^. Restated in a more sophisticated way, for a given 



loop 17 the map $^ reconstructs the Bloch bundle £^{r}) as the pullback [Hus66| of the 
reference vector bundle S^im, i. e. £5(17) ^ ^*(.£s^'«)- 



T)^ 



4.3 The case of two internal degrees of freedom 

Let us now specialize to the simplest, but still non-trivial case of a system with only two 
internal degrees of freedom, i. e. m = 1. In this case we identify the T,i,T,2,T,2 with the 
Pauh matrices CTi, ctj, CT3. With this choice the reference projector ( |4.5P reads 



r. rr. ^ 1 / 1 - COS -e"'^ siu „ , ,_„, 

Ps^(^,^)=2 _e+..,i,0 l + ,„30 I, (4.10) 



and the inverse transforms ( |4.9| l can be equivalently written as 



'h„olik,t)\ 
^''^^'^) = ^^^^°M^^J' W.ll) 

'^2(^,0 
(/3^(fc, t) = arctan ' 



It is possible to provide an explicit local trivialization for the complex line bundle ^§2: 
One needs at least two charts to cover §^, and we shall always choose what we call an 
AfS covering: let Cj^ and Cg be two closed, contractible and mutually disjoint sets with 
non-empty interior which contain only either the north pole {6 = 0) or the south pole 
{6 = Ti) pole (denoted by U and S, respectively). Then we define Uj^ := S^ \ Cg and 
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4.3 The case of two internal degrees of freedom 



Us := §^ \ Cj^ which serve as our open covering of the sphere. From the definition it 
follows that Uj^ f\Ug^%. On the northern hemisphere, we may choose 



-IV 



sm; 



■cos 2 



*^(0,^):= _e2 (4.12) 



as local section, i. e. {0,ip) -^ "^j^^O,^) is continuous and span ranPg2(0,c^) for all 
(6, (/?) e Uj^. However, this parametrization cannot be extended unambiguously up to the 
south pole S since *jy-(5) = (e~"^, O) depends non-trivially on ^p. But the parametrization 



sin — 



is well-defined at S (but not at W), and (0, (/;) -^ *5(0, i/') is a local section on Ug. These 
local sections are glued together with the transition function 



We recall that all complex line bundles f ^ X are completely classified by the first Chem 
class Ci(f ) e H^[X,Z). Moreover, the Chern-Weil theory allows us to consider c^^S) as a 
differential 2-form (modulo exact forms), i. e. Ci(£) e H^j^CX) [MS74||. In the particular 
case of the reference line bundle ti : £§2 — > S^ we can provide an explicit representative 



for Ci(£g2) simply following [BT82 pp. 71-75] which uses the identification between top 
Chem classes of complex vector bundles and Euler classes of the corresponding realifica- 
tions. This gives 

Ci(fg2) = dC 



where ( := {^j^, ^5} is the collection of two local 1-forms defined by 



~i2^^'sdlngAA5 on U^nUg 

r .= / + I^(^~^'s)dlngA^5 on U^nU^ 

^' [0 on Cs 

where Xs i^ ^ny smooth function supported in U^ . The first Chern number of the line 
bundle £^§2 is by definition 



C(fs2):= Ci(fs2)= dC. 

Now, let /Cjv- ^ Cjv" an open set such that dlCj^ is a sufficiently regular closed path around 
Af. Moreover, we assume without loss of generality that XsldKj^ ~ ■*■■ Then we can use 
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4 Two-band systems 



the Stoke's theorem to write: 



JKj^ js^\Kj^ J J J 

dKf^ dKf^ dKj^ 

Using the explicit form of the local forms i^^^^ and i^g and of the transition function gj^g 
we arrive at the important formula 

C(£s2)=^i>d(^ (4.13) 

A 

where A is any closed regular path which goes around the north pole M in the counter- 
clockwise direction (i. e. the direction of increasing t^). Since y> e [0, 27r] wraps around 
once, we obtain C(£g2) = 1, i. e. £§2 is non-trivial. 

Now we can return to the study of the Bloch bundle S^{rj) c^ $*(£g2). The topology 
of this line bundle is fully described by the first Chern class Cj (£"3(17)) which can be 
reconstructed from CjCfgi) with the map $*. In fact, t\ye functoriality (or naturality) of 
the characteristic classes [ |MS74P implies 

ci (£3(7,)) = ci ($;(£s2)) = $; ci(£s2) (4.14) 

where the last $* denotes the induced group morphism between the cohomology groups 
h2(§2^Z) and H^(J'^+\Z') (resp. between the de Rham groups H^^^CS^) and H^j^CT'^+i)). 
The pullback structure of the Bloch line bundle £b(^) allows us to define a local trivial- 
ization inherited from £§2. The relevant covering of the base space T'*"'"^ is given by 

^AA,, := ^-\i^^), Us,^ := %\Us). 

Similarly, we will define the north and south pole variety as 

N^ := $-i({AA}) = {(fc, t) e T'^+i | 0^(/c, t) = O}, 
5, := 4.-i({5}) = {{k, t) e T'^+i | 0^(/c, t) = ti}. 

On U^.^ and ZY^ ,,, the compositions *(j,^(fc, t) := ^j o $^(fc, t) = *j(0^(fc, t), (/3^(fc, t)), 
with tt = A/", 5 respectively, are local frames which are glued together through the tran- 
sition function gTV's.r; '■= Sms ° *77 = e"*""^". This give us an immediate criterium for the 
triviality of the Bloch bundle £b(-V)'- 

Theorem 4.2 //7V^ = or 5^ = 0, then the Bloch line bundle S^ir]) is trivial. 

So let us proceed to compute ( |4.14| ): We recall that the T^^ '-^ T''^^ with 1 ^ j < n ^ 



d + 1, are the generators of the homology group HjCT''"'"^) described in Section 3.3 The 
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integration of the 2-form Cj (i?b(^)) over these two-dimensional surfaces produces a set 
of independent Chern numbers 

C;,„([r?]):= f Ci(fB(r,))=f ^;c,[£^2l 

The Chern class Ci[£^[r])) can be realized as the exterior derivative of the locally defined 
1-forms {d(j^jj,d(g.^} where d^j^^ := $*(d(^j) are defined onZ^j^j with ^ = Af,S. 

Now we can pick an open set ICj^,^ Z) Cj^^j := ^~^(_Cj^) which contains the north pole 
variety A/"^ with a sufficiently regular boundary dKj^.^ c Uj^^^ PiUg.^. If we can embed 
'^^n ^ '^'''^^ su'^h that Tj^ n ICj^_^ = 0, then Stoke's theorem implies Cj_n([r]]) = 0. Put 
another way, the Chern class Ci(^£^[r])) behaves like the exact 2-form dij^^ on T^^. This 
gives yet another simple criterion: if 7V^ n T^^ = or 5^ n T^^ = 0, then C^ ,,([17]) = 0. 

However, if both the intersections are non-empty, then we exploit 

d (T^^ n /C^, J = d [Tl^ \ (T^^ n /C^,^)) = t^^ n 5/C^,^ 
and dlngjV5 ° $,, = idi/?^ to obtain a description of C^ ^([tj]) as a winding number, 

C;,„([rj])= A (Cs.r, - Ca^,,,) = ^ i d(^^. (4.15) 

9(Tj?„nJC^,„) Ai,„ 

Here A^ „ is any closed regular path which goes around the north singular set Af^ n T^ ^ in 
the counter-clockwise direction. 



5 The uniaxial strain model 

We now compute ( |4.15| ) for the simplest kind of model for piezo effects in graphene, 
the uniaxial strain model. It allows for nearest neighbor interactions with three hopping 
parameters qg, qi and qj as well as a stagger parameter q^: 

V^qo -If^fri + qiSj +q2S2 -q3-lf2(r) J 

In principle, the parameters can be complex, but to ensure selfadjointness of H, q^ needs 
to be real. 

To understand the action of this hamiltonian, we recall Figure [T| the honeycomb lattice 
consists of two sublattices, and each (white/black) atom has three (black/wite) nearest 
neighbors. The offdiagonal terms describe nearest-neighbor hopping. If an electron is 



31 



5 The uniaxial strain model 



initially located at the white atom at y, it hops to its neighboring black sites located at y, 
y — Yi and 7 — 72 with rates qo, qi and qj (green vectors) . Applying stagger qa 7^ implies 
that white and black atoms within a unit cell are no longer equivalent, thus breaking the 
inversion symmetry. 

5.1 Description of the gapped parameter space 

Instead of the periodic two-band operator ( |5.ip , we will study its counterpart in the Bril- 
louin algebra Ap^,-, namely 

H(fc; q) = Re (c7(fc; q)) a, + Im (c7(fc; q)) 0-2 + qs CT3, (5.2) 

where we have introduced the function 

Cj(fc;q) = qo + qie-''=i+q2e-*^ 

and the short-hand notation q = (qo,qi,q2j'J3) e C^ x R, fc = (fcj,fc2) ^ T^. For fixed 
(fc;q), the eigenvalues of H(fc;q) are £±(A:;q) = ±^/q^+\m(qyk^, and hence the total 
spectrum 

CT(H(q)) := y C7(H(fc;q)) =/_(q)U/+(q) 

consists of two bands /±(q) := ran£-|-(-;q). Clearly cr(H(q)) is symmetric around E = 
and has a gap there if and only if £+(fc;q) > holds for all fc e T^. This allows us to define 
the gapped parameter space (C'^ x K)q = (C^ x K) \ Q^g at Fermi energy E = where 



Proposition 5.1 ( flHKN-l-Oe] ) 



Qng:={(0,qi,q2.0) | |qil = |q2l}u 

u{(qo,qi,q2,0) I qo#0, ||9i/qo|-|«2/9ol| ^l^l^V'iol + l'iV'iolll- 
Proof The proof is straightforward: from the form of the eigenvalues, £±(fc; q), it is clear 
that qg = is a necessary condition for there to be no gap. So let q^ = 0. Then |c7(q; fc)| = 
is equivalent to uj[q;k) = 0. In case qg = 0, then Iq^l = Iqjl is a necessary, and 
indeed sufficient condition (the points of intersection depend on the phase of the complex 
numbers qj and q2). In the remaining case, qo 7^ can be factored out and we can rewrite 
C7(q;fc) = as 

90 % 
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5.2 The Chern numbers 



and taking the square of the modulus of both sides yields 

2 



2 

+ 



+ 2 Re I ^e+'^^.-<^-^ 



The right-hand side takes values in [(pj — P2)^,(Pi +P2)^]) Pj •= hVioL and thus the 
above equation has a solution if and only if q is in the second region of Q^g. □ 

The gapped parameter space is independent of phase and sign of the hopping parameters. 
Moreover, as the model is trivial if two of the three hopping parameters vanish, we can 
set qo = 1 by choosing a suitable energy scale for the hamiltonian ( |5.2| ). Hence, without 
loss of generality, we can restrict our attention to the parameter space 

Q = {1} X [0, 2] X [0, 2] X [-1, 1] c C^ X K. (5.3) 

When necessary, we will smoothen the corners of this cube so that Q is a C°° -manifold. 
(But even with sharp corners 7iJ(Q) ^ ttiCQ) holds which is all we really need.) We have 
sketched the parameter space of gapped configurations Qq = Q\Qng in Figure[2|(b). Then 
one immediately deduces the following result: 

Proposition 5.2 (Fundamental group of Qq) The fundamental group 7ii(Qo) ^Z'^is gen- 
erated by 

Tji(t)= (l-l-e cost,0, -e sint), rjjCt) = (0, 1 -I- e cos t,-e sin t), (5.4) 

for some e small enough. 

A detailed proof of this quite obvious fact would be unnecessarily technical, so let us just 
sketch the basic idea using Figure [2] (b) : It is clear that 17 ^ and 172 each generate distinct 
fundamental loops. One may think that loops which wrap around the gapless region and 
run from qj = (3/2, 1/4, 0) to q^ = (1/4, V2, 0) and back represent yet another distinct class 
of loops. However, these loops can be deformed into a composition of loops equivalent to 
172-171- 

5.2 The Chern numbers 

On Q, the three prefactors ( |5.2| | of the Pauli matrices simplify to 



/ii(fc,q) := Re (tD-(q;fc)) = 1 -t-q^ cosfcj -|-q2 cosfc2 

?i2(fc;q) := Im(c7(q;fc)) = qj^ sinfcj -|-q2 sinfc2 (5.5) 

h3(fc;q) = q3 
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5 The uniaxial strain model 



The symmetry of the model manifests itself in 



5.2 



Hence, 



where the rj^ are the generators of the fundamental groups from Proposition , 

it suffices to calculate the Chern numbers for one of the generators. The main purpose of 

this section is to prove the following 

Theorem 5.3 Let 17 e C(S^,Qq) be a loop and [17] = n^ [tjj] + Hj [172] e tzi^Qq) its equiv- 
alence class where the r]j are given by ( |5.4| ). Then the matrix of Chern numbers (cf. Sec- 
tion\. 



3.31 is given by 



C([rj]) = niC([rji]) + n2C([rj2]) 

f o|iA f 



I 
V -f ~o"; b 



+ n, 







We follow the ideas put forth in Section 4.3 in order to compute Ci([i7i]) component-by- 



component. In view of ( |5.5| ) and ( |4.1lD , the angle coordinates are 

I —e sin t 

On^(.ki,p^,t) = aTccos 



\/l + 2(l + e cost) cosfcj-l-Cl-l-e cost)^-!-;?^ 



sin^ t 



f (1 -l-e cost) sinfcj 

(^n (fci , K, = arctan ; ; — 

^''^ ^-^^ Vl-l-(l-l-ecost)cosfci 

The north and south pole varieties are determined by 6^ : for A/^ and S^ , we need to 
solve 



e sin t = =F V 1 -I- 2(1 -I- e cos t) cos fci -I- (1 -I- e cos t)^ -I- e^ sin^ t, 
and a short computation yields 

AA^^ = {(-7r,fc2,-f) I k2&[-7i,+n)}, 
5^^ = {(-7r,fc2, + f) I k2&[-7i,+n)}. 
Our first task is to verify whether the intersections of the tori 

T?, = {(fci,fc2,0eT3 I k2 = *&[-n,+Ti)} 



T: 



T' 



1,3 

i2 
2,3 

i2 



{(fci,fc2,t)eT3 I fci=*e[-7r,+7r)} 
{(fci,fc2,t)eT3 I t = *&[-7i,+n)} 
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5.2 The Chern numbers 



with the north and south pole varieties are empty, because if at least one of them is, the 
Chern number has to be 0. The intersections with the first torus consists only of a single 
point, 

Af,,ml, = {(-^,*,-|)}, s^^ml, = {{-n,*, f )} . 

Whether the intersections with T^ 3 are empty or A/"^ and S^ themselves depends on the 
choice of embedding point: 

From a topological perspective, the embedding point (i. e. the choice of fcj = *) is not 
important, so we conclude C23(_[r]i]) = 0. For the sake of completeness, we will postpone 
the explicit computation for the case * = tt to the end of the section. In case of the last 
torus, at least one of the intersections is always empty independently of the choice of *, 
and thus Q 2([i7i]) = 0. 

To compute Ci3([tji]), we pick an open set JCj^^^^ which contains pj^ = (— ^,*, — f ) 
so that the boundary dlCj^.^^ is parametrized by the circle 

71,3(5) = pj^ + 5 (coss, 0, sins) 

for some 5 > small enough. The exterior derivative 

depends on 

(1 + e cos t)(l + e cos t + cos k^) 

^Jc, Vri, (fcl> = ~T , nr-j , r^ ; , r-j , T^' 

' " 1 + 2(1 + e cost) cosfcj + (1 + e cos ty 

e sinfcj sint 
t^r,S 1' '^-' ~ 1 + 2 (1 + e COS t) cosfcj + (1 + e cos t)^ ' 

Expanding the partial derivatives evaluated at 7i_3(s) for small 5, we obtain 

—5e sins + 0(5^) 1 e sins 

3k, V„i (roi(5)) = - TTr 2 _L 2 ■ 2 ^^,nrK3^ = ~ ^ 2 _l 2 ■ 2 + ^(1)' 

o^(cos^s + e^ sm sj +0(5'') ocos^s + e^sms 

—5e coss + 0(5^) 1 e coss 

' 5^(cos^s + £:^ sm^s) +0(5-^) o cos^ s + e^ sm^ s 
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5 The uniaxial strain model 



Now we plug in the parametrization into equation ( |4.15| ) to obtain 

Ci.3([r}i])=— <j) ^'^v, = ^\ dsVi/j^^(ri_3(s))-|-i3(s) 






2ji , 

as 
27T jg cos^s + e^ sin^s 



To arrive at the last equality, we note that the leading-order integral can be computed 
explicitly and that Chern numbers need to be integers. 

For the benefit of the reader, we will double-check that C2_3([i7i]) = exphcitly, even if 
we choose the "bad" embedding point * = —n. Here, the singular string A/"^ c T^ 3 is fully 
contained in the torus. Pick an open set ICj^ ^^ Z) A/^^ so that 3/C^ ^^ can be parametrized 
by 

r2Vs):=(-7r,s,-f±5) 

for 5 > small enough. Then from 72 3^^^) ~ ^^' ^'^^ '^^'^ ^h'^vi ifis^^^) ~ ^' ^^ ^^^' 
elude that 

which implies that C2_3([rji]) vanishes. 

5.3 The effect of perturbations 

We are interested in discussing the effect of a small perturbation XV (periodic or random) 



to the operator ( |5.2p . As in Section 1.3 we will assume A <: 1 and ||V||^ = 1. 

The assumption that there exists a global gap of size at least g for the hamiltonian ( |5.2| ) 
is equivalent to 

ql + \co[q;k)\^>^, VfceT^. 

This condition is of course realized if |q3 1 > |. For (qi,q2,0) e Qg, the above condition is 
realized if one of the following inequalities is verified 

S S S 

q2<qi-i--, q2>qi + i+2' q2<-qi + i-2- 

The space Qq^ c Qq obtained by the configurations which verify the above relations is rep- 
resented in Figure Is] It is clear that Qq^ -^ Qq when g ^ 0, i. e. Qq^ deformation retracts 



to the thin set Qq in a topological sense [ Hat02 1 . This means that Qq „ is homotopically 



equivalent to Qg (at least for g < 2) and so 7ii(Qo ,,) = tzi^Qo)- 
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Any Tj e Qog can also be considered as a loop in Qq and its homotopy class can be 
expressed as 17 = ni[i7i] + n2[r]2]- Moreover, for all A e [0, A,] with A^ < | we can 
compute the Chern class C^>,a( [''?]) for the perturbed hamiltonian Hx^^iM- The homotopy 
equivalence of the spectral projections (Proposition 2.41 implies Cj„;i([i7]) = C^_„([i7]) 
agrees with the Chern numbers from the periodic case. 



A Symmetric classes for two-band periodic systems 

Let us consider a two-band periodic operator in the algebra Qlp^j. = 6 ® MatjmCC) of the 
form 



H 



2m+l 



(A.1) 



where f)^ = [)* e ©. In definition ( |A.1| ) there is no loss of generality compared to ( |4.1 
since H and H — [}o ® 12m have the same spectral projections. 
An explicit realization for the Clifford matrices E^ is given by 



CTi 



> C7, 



iCTl 



i a-. 



) a^ 



i a-. 



) (To 



-■2771 



±2'^ -12'^ ±2 
1, (8) 1,® 1 



ICTi 



•■2 

0"3 



2 W ±2 ^ 
C73 18) 0'3 I 



i (To 



(A.2) 



where each term is a tensor products of m Pauli matrices Ui, CT2, 0-3. This choice is essen- 
tially unique, up to unitary equivalences [Lee48] . Since a^, a^ are real matrices and o"2 is 
purely imaginary it follows that T,j = (— ly^^ T,j. 

The peculiarity of these models resides in the existence of a symmetry which is reflected 
in the spectral properties. First of all, let us recall the complex conjugation which acts 
antilinearly as C : * -^ * for * e i^(X) ® C' (r is arbitrary). Moreover, we will also need 
the parity operator p on ^^(F) which in particular acts by conjugation on selfdajoint (j e © 
aspt)p = C[)C. Moreover p = p* and p^ = l|2(-pj. 

Let us start with the odd case m = 2v — 1 and introduce the matrix 



:— S]^S3 . . .S4,,_i, 



e* = e- 



e^ 



K-1712- 



(A.3) 
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Since 9 is a product of 2v matrices (with odd indices), exploiting the CHfford algebra 
relations one obtains &T,j = [—iyT,j&. If one introduces the unitary operator Uq := 
piS>Q, with involution property U^ = (—1)^' 1-^, one can verify that 

UqHUq = -CHC. (A. 4) 

The relation ( |A.4| ) says that the operator H has a particle-hole (PH) symmetry in the 



language of [ SRF+08 1 



For the even case m = 2v one introduces the matrix 



T := es. 



'4v+l) 






T^ = (-1)" l2,„ 



(A.5) 



Since T is a product of 2v+l matrices (with odd indices), one obtains T S^ = [—ly'^^T.jT. 
The unitary operator t/y := p T, with involution property U^ = (—1)' 1^, provides 



UrHU* 



CHC. 



.^-^..^. (A.6) 

The relation ( 1A.6D says that the operator H has a time-reversal (TR) symmetry again ac- 



cording to [ SRF+08 1 



One has the following table of symmetries for H depending on the spinorial dimension: 



m 


V 


PH 


TR 


AZ 


even 


even 





+ 1 


AI 


even 


odd 





-1 


All 


odd 


even 


+1 





D 


odd 


odd 


-1 





C 



Table 1 : Table of the discrete symmetries for two-band periodic hamiltonian H of type 
dA.lD . The classification depend on the dimension of the internal degrees of 
freedom through the parity of m and v := [(m + l)/2j. The sign of the PH- 
symmetry (resp. TR-symmetry) is the sign of U^ (resp. U^). The means absence 
of symmetry. The colum n AZ provides th e Cartan's label according to the Altland- 
Zirnbauer classification []AZ97[ [sRF+08 | . 
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